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ABSTRACT 


A  linearized  perturbation  analysis  is  used  to  determine  if  insta¬ 
bilities  can  occur  in  MHD  accelerators.  When  the  electrical  pertur¬ 
bations  are  assumed  unaffected  by  the  presence  of  conducting  walls, 
the  enthalpy  wave  has  the  largest  growth  rate.  The  growth  rate  de¬ 
pends  upon  the  orientation  of  the  current  vector  to  the  wave  front  and 
the  location  within  the  accelerator.  When  the  electrical  perturbations 
are  influenced  by  the  presence  of  the  slanted  conducting  walls,  the 
upstream-moving  magneto -acoustic  wave  is  amplified.  When  steady- 
state  gradients  are  neglected,  the  growth  rate  changes  only  slightly  in 
both  cases.  An  elevated  electron  temperature  increases  the  enthalpy 
wave  growth  rate  in  the  absence  of  conducting  walls,  but  has  little 
effect  when  the  conducting  walls  control  the  electrical  perturbations. 
The  linearized  analysis  indicates  that  instabilities  can  exist  in  MHD 
accelerators  and  the  growth  rate  may  be  large  in  some  cases. 
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NOMENCLATURE 


A 

Area 

Al,  A-2 

Defined  after  Eq.  (51)  for  one -dimensional 
perturbations 

Ai,  A2,  A3 

Defined  in  Eq.  (86)  for  two-dimensional 
perturbations 

a 

Constant  defined  in  Eq.  (28) 

a0 

Speed  of  sound 

al» a2»  a3 

Defined  in -Eq.  (1-9) 

B 

Magnetic  flux  density 

Bi 

Defined  in  Eq.  (28) 
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Blj  B2,  B3 


b 

bl>  b2j  b3 
bj  through  b6 
Cj  through  Cg 
CAi  through  CA4 
CB^  through  CB4 
D 
d 
E 

El,E2,E3 

e 

F 

FXJ  Fy 

Gi,G2,G3 

H 

h 

I 

J 

3 

k 

L 

M 


m 

N,Ni,N2 

Nef 

n 

P 

P 


Defined  after  Eq.  (51)  for  one -dimensional  case 
and  in  Eq.  (86)  for  the  two-dimensional  case 

Defined  in  Eq.  (28) 

Defined  in  Eq.  (1-9)  for  one -dimensional  case 
Defined  in  Eq.  (86)  for  the  two-dimensional  case 
Coefficients  in  Dispersion  Relation  Eq.  (52) 
Coefficients  in  Dispersion  Relation  Eq.  (87) 
Coefficients  in  Dispersion  Relation  Eq.  (87) 
Defined  in  Eq.  (62) 

Channel  height 
Electric  field 
Defined  in  Eq.  (86) 

Perturbed  electric  field 
Defined  in  Eq.  (61) 

Defined  in  Eq.  (86) 

Defined  after  Eq.  (51) 

Total  enthalpy 

Static  enthalpy 

Total  current 

Current  density 

Perturbed  current  density 

Wave  number  or  Boltzmann's  constant 

Characteristic  length  in  x  direction 

Mach  number 

Particle  mass 

Defined  after  Eq.  (51) 

Defined  in  Eq.  (86) 

Particle  number  density 
Defined  in  Eq.  (86) 

Pressure 


Vll 
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q2  through  q4 
R 


Sm 

T 

V 

Vi 

u,v 


w 

x,y,  z 

a 

& 

@ef 

7 

7i  through  73 

'  5S 

6eff 

€ 

6 

\ 

v 

Si* §2i §3 
P 

CT 

CTef 


Defined  after  Eq.  (1-4)  for  Case  I  and  after 
Eq.  (56)  for  Case  II 

Defined  after  Eq.  (45)  for  Case  I 

Gas  constant 

Defined  in  Eq.  (60) 

Seed  fraction  by  mole 

Temperature 

Perturbed  electric  potential 
Seed  ionization  energy 

Velocity  components  in  x,  y  directions,  respec¬ 
tively 

Enthalpy  exponent  defined  after  Eq.  (66) 
Cartesian  coordinates 
Ionization  fraction 
Hall  coefficient 

Reduced  Hall  coefficient  =  fiQ/(l  +  /302) 

Gas  specific  heat  ratio 

Defined  after  Eq.  (46)  for  Case  I  and  after 
Eq.  (56)  for  Case  H 

Energy  loss  factor  of  electrons  interacting  with 
species  s 

Effective  energy  loss  factor  of  the  gas 
Angle  between  current  vector  and  wave  front 
Wall  slant  angle 
Disturbance  wavelength 
Defined  after  Eq.  (1-6) 

Collision  frequency 
Defined  after  Eq.  (1-6) 

Density 

Electrical  conductivity 

Reduced  electrical  conductivity  =  aQ/(l  +  |30^) 


viii 


A  EDC-TR -70-290 


r  Average  time  between  collisions 

®  Tangent  of  wall  slant  angle 

<l>  Defined  after  Eq.  (1-4) 

< i i  Perturbed  velocity  potential  function 

tjj  Perturbed  current  density  potential  function 

10  Complex  disturbance  frequency  =  +  ito2 

SUBSCRIPTS 

1 
e 
f 
i 

K 
s 

X 

y 

o 


Perturbation  quantity 

Electrons 

Final  value 

Initial  value 

Potassium 

Species  other  than  electrons 
Component  in  x  direction 
Component  in  y  direction 
Steady-state  quantity 
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SECTION  I 
INTRODUCTION 


The  prediction  of  magnetohydrodynamic  (MHD)  accelerator  perform¬ 
ance  requires  an  understanding  of  the  basic  flow  phenomena  and  param¬ 
eters  which  affect  their  operation.  The  presence  of  any  disturbance  which 
is  significantly  amplified  can  greatly  affect  the  operation  of  these  devices. 
Previous  investigations  considering  the  growth  of  disturbances  in  MHD 
generators  have  indicated  that  the  growth  can  be  significant  in  some 
cases.  The  purpose  of  this  investigation  is  to  determine  if  disturbances 
can  grow  in  MHD  accelerators  resulting  in  an  instability  which  can  affect 
the  performance  of  the  device. 

There  are  three  types  of  disturbances  which  may  result  in  an  insta¬ 
bility  in  the  high  density,  highly  resistive  plasmas  utilized  in  MHD  accel¬ 
erators.  The  "electrothermal"  (ionization)  instability,  first  postulated 
by  Kerrebrock  (Ref.  1),  is  the  result  of  preferential  heating  of  the  elec¬ 
tron  gas  with  the  heavy  gas  species  unaffected.  One  necessary  criterion 
for  the  amplification  of  these  "electrothermal  waves"  is  that  the  electron 
gas  be  loosely  coupled  to  the  heavy  gas  (Refs.  1  and  2).  The  coupling 
between  the  electrons  and  heavy  species  in  MHD  accelerators  using 
diatomic  gases  is  quite  strong  and,  hence,  the  electrothermal  disturbance 
is  not  expected  to  lead  to  instabilities. 

The  enthalpy  wave  is  similar  to  the  electrothermal  wave  but  involves 
localized  heating  of  the  entire  gas.  A  localized  "hot  spot"  in  the  gas  is 
heated  further  because  of  the  local  increase  in  conductivity  and  hence  in¬ 
creased  joule  heating.  The  growth  rate  of  an  enthalpy  wave  will  be  much 
less  than  that  of  an  electrothermal  wave  because  of  the  larger  specific 
heat  of  the  gas.  However,  in  contrast  to  the  electrothermal  wave,  the 
enthalpy  wave  can  exist  in  the  MHD  accelerators  of  interest. 

The  amplification  of  the  magneto -acoustic  disturbance  is  directly 
dependent  on  the  fluctuations  of  the  thermodynamic  properties  of  the  gas 
(Ref.  3)  and  has  a  large  growth  when  the  electrons  and  heavy  species  are 
strongly  coupled  (Ref.  2).  This  type  of  instability,  first  postulated  by 
Velikhov  (Ref.  4),  depends  directly  on  the  Hall  effect  in  the  first  approxi¬ 
mation.  A  simple  physical  model  of  this  disturbance  can  be  constructed 
in  the  following  manner  (Refs.  2  and  3).  Consider  the  case  of  the  total 
steady-state  current  in  the  y  direction  and  the  magnetic  field,  B,  in  the 
z  direction  as  shown  in  Fig.  1. 
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In  the  presence  of  a  positive  density  perturbation  propagating 
parallel  or  antiparallel  toTp,  an  induced  current,  0-o(^  x  5),  will  flow. 
The  resulting  force  a0(u  xB)  will  always  oppose  the  propagation 
and,  hence,  damp  the  disturbance.  However,  the  positive  density  per¬ 
turbation  results  in  a  local  decrease  in  the  Hall  parameter,  /3.  The 
decrease  in  d  means  that  the  x  component  of  the  electric  field  has  too 
large  a  value  for  the  local  conditions  and  therefore  a  perturbation 
current,  J^,  will  flow  in  the  positive  x  direction.  This  current  inter¬ 
acting  with  the  magnetic  field  will  amplify  the  disturbance  moving  anti¬ 
parallel  to  JQ  and  attenuate  the  disturbance  moving  parallel  to  JQ.  The 
condition  for  instability,  given  initially  by  Velikhov  (Ref.  4),  is  that 
Jh  >  CTo  u  B  when  there  are  fluctuations  only  in 

Detailed  analyses  (Ref.  3),  including  perturbations  in  conductivity, 
pressure,  density,  velocity,  and  Hall  parameter,  indicate  that,  de¬ 
pending  on  the  frequency,  the  initial  steady  state  of  the  plasma,  and 
the  orientation  of  the  wave  vector,  waves  can  be  amplified  at  a  rate 
significantly  higher  than  that  indicated  by  Velikhov.  The  analyses 
(Refs.  5,  6,  and  7)  of  the  transverse  and  longitudinal  waves  in  a  gener¬ 
ator  operating  in  the  Faraday  mode  indicate  that  the  growth  rate  is 
small  and  the  instability  can  be  controlled  by  external  circuitry.  How¬ 
ever,  the  growth  of  axial  waves  in  a  Hall -type  generator  was  found  to 
be  significant  (Ref.  8).  The  growth  rates  of  the  waves  have,  in  general, 
been  found  to  increase  with  decreasing  pressure  (Ref.  9)  or  increasing 
specific  heat  ratio  (Refs.  8  and  10).  The  results  of  these  investigations 
indicate  that  instabilities  can  exist  in  MHD  generators. 

On  the  basis  of  the  substantial  evidence  that  instabilities  can  exist 
in  MHD  generators,  there  is  reason  to  believe  that  instabilities  can 
also  exist  in  accelerators.  The  investigation  of  magneto-acoustic  and 
enthalpy  disturbances  in  MHD  accelerators  is  the  primary  concern  of 
this  report.  The  method  of  analysis  will  be  to  assume  that  the'  steady- 
state  flow  is  given  by  the  solution  of  the  usual  quasi- one -dimensional 
gas  dynamic  equations  with  MHD  terms.  The  flow  is  subjected  to  a 
small  perturbation,  and  the  behavior  of  this  distrubance  is  investigated 
using  a  linearized  analysis. 


SECTION  II 

MATHEMATICAL  MODEL 


The  equations  describing  the  flow  of  an  inviscid,  nonheat¬ 
conducting,  slightly  ionized  gas  in  the  presence  of  electric  and  mag¬ 
netic  fields  can  be  written  in  the  form  (Ref.  11): 
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Mass : 


dp  — 

— ^  +  7  •  pv  =  0 

at 


(i) 


“•  “• 

Momentum:  p  —  +  pv  •  vv 

dt 


-  vp  +  J  x  B 


(2) 


_  dH  —  dp  “T 

Energy:  p  —  +  pv  *  VH - —  =  J  -  E 

dt  dt 


(3) 


■where  p  is  the  density  .^v  is  the  velocity,  p  is  the  pressure,  J  is  the 
total  current  density,  B  is  the  applied  magnetic  field,  E  is  the  sta¬ 
tionary  electric  field,  and  H  is  the  total  enthalpy  given  by 


H  =  h  + 


v  « 


(4) 


where  h  is  the  static  enthalpy.  The  electric  current  density  is  given 
by  Ohm's  law,  including  the  Hall  effect,  as 


1  +  0 


2 


j  -k 

E  +  v*B 


-»( 


ExB  +  (vxB) 


xS]J 


(5) 


where  cr  is  the  electrical  conductivity  and  J3  is  the  Hall  coefficient. 


2.1  STEADY-STATE  SOLUTION 


The  steady  accelerator  flow  is  assumed  to  be  given  by  the  quasi- 
one -dimensional  gas  dynamic  equations  including  MHD  effects.  Using 
the  quasi-one-dimensional  approximations  in  Eqs.  (1),  (2),  and  (3), 
the  equations  for  the  steady  state  are 


p  u  A 
yo  o 

• 

=  m 

(6) 

dx 

dpo 

■ - °  +  J  B 

dx  y 

(7) 

dH 

_ o 

=  J  E  +  J  E 

(8) 

dx 

xx  y  y 
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where  u  is  the  velocity  in  the  x  direction  and  m  is  the  mass  flow  rate. 

An  MHD  accelerator  channel  having  slanted  electrodes  and  con¬ 
nected  in  the  two-terminal  mode  (Fig.  2  and  Ref.  12)  is  to  be  considered. 

METALLIC 


Fig.  2  Slant-Wall  Accelerator  Connected  for  Two-Terminal  Operation 


For  this  case  the  generalized  Ohm's  law  can  be  written  in  component 
form  as 

j - 12 — _  /E  -  S  (E  -  u  B)}  (9 

1  +  P0  L  7  J 


5 


AEDC-TR -70-290 


where  4>0  =  tan  0,  I  is  the  total  current  supplied  to  the  channel,  A  is 
the  channel  cross-sectional  area,  and  6  is  the  channel  slant-wall  angle. 
The  voltage  drop  at  the  electrode  surface  has  been  neglected  in  this 
analysis.  Using  Eqs.  (11)  and  (12)  in  Eqs.  (9)  and  (10)  results  in  the 
following  expressions  for  Jx,  Jy,  and  Ex: 


i  ‘  Mo, 

*  (T7TT)  + 


f 

- 2 -  o  u  B 

1  +  t  a  °  ° 


(13) 


T  i  +  0  a  u  B 

±  r  Q  Ox  o  o 
A  \  .  2}  "  i  +  4  2 


1  +  * 


(14) 


E. 


(I/fl)(l  +  30a)  -  °oaoB(3o  -  l0) 


°o(1  +  *o2) 


(15) 


Solutions  to  Eqs.  (6)  through  (8)  and  (13)  through  (15)  for  given  values 
of  I,  A,  and  B,  which  may  vary  with  x,  and  with  specified  initial  con¬ 
ditions  have  been  obtained.  The  solution  shown  in  Fig.  3  for  potassium- 
seeded  air  in  chemical  equilibrium  was  used  for  the  stability  analysis 
presented  here.  Other  operating  conditions  were  considered,  but  the 
general  criteria  for  instability  did  not  change  from  the  results  obtained 
utilizing  this  solution. 


2.2  PERTURBATION  EQUATIONS 


The  response  of  the  plasma  to  small  disturbances  will  be  investi¬ 
gated  using  a  small  perturbation  analysis.  The  externally  applied 
magnetic  field  is  assumed  constant  with  respect  to  the  perturbation 
and  induced  magnetic  fields  are  neglected  since  the  magnetic  Reynolds 
number  is  small.  The  velocity,  thermodynamic,  and  electrical  param¬ 
eters  are  assumed  to  undergo  a  small  perturbation;  that  is 


V  =  VQ  +  Vi  =  i(uQ  +  Ui)  +  j  (v  i) 

7  =  7o  +T  -  ^(jx  +  jx>  +  ?(jy  +  jy) 

E  =  Eq  +  e  =  i(Ex  +  ex)  +  j (Ey  +  ey) 


(16) 
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Fig.  3  Steady-State  Solution  for  a  Slant-Wall  Accelerator 
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Substituting  the  perturbed  variables  into  Eqs.  (1),  (2),  (3),  (9),  and 
(10)  and  neglecting  higher-order  terms  results  in  the  following  per¬ 
turbation  equations: 


Continuity: 

^Pi  ^Pi  duQ  SpQ  Sui  5vi 

sr  +  “o  ar +  P1  ar +  U1  sr +  po(3r +  it5  " 0  <17> 

x- Momentum: 

3ui  ^ui  Sun  pi  ui  Spi 

p  -  +  p  u  +  p  u  — —  ( —  +  — )  =  -  — —  +  j  B  (18) 

°  3t  °  °  Sx  °  °  Sx  PQ  Uo  ^X  ^ 

y  -  Momentum : 

Svi  5vx  ^Pi 

po  St  Po  o  Sx  Jx 

Energy: 


(19) 


dHi  Sho  ui  pi  ^Pi 

p  -  - —  +  p  u  - —  +  p  ii  - —  ( —  +  — )  -  

°  St  °  °  Sx  °  °  Sx  uQ  p0  St 


(20) 


E0-j  +  e.JQ 


Ohm's  Law: 

Jx  -  aef  1  ex  +  V 


•I  e  +  ViB  -  p  (e  -  uiB)  -  Pi(E  -  u  B)  f 
\  x  *  Ko  y  M  y  o  J 

+  Jx-[a^  "  2  gefPlJ> 


(21) 


jy  =  CTef  {ey  "  UlB  +  Po(ex  +  VlB)  +  3lEx} 


+  J  \  —  —  2  3  -3  i  !■ 

U^o  ef  J 


where 


(22) 


ef 


1  +  P, 


ef 


1  +  P, 
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This  general  form  of  the  perturbed  Ohm's  law  will  be  utilized  in  later 
sections  when  special  forms  of  the  perturbed  equations  are  considered. 

The  four  partial  differential  equations  (Eqs.  (17)  through  (20)).  and 
the  two  algebraic  equations  (Eqs.  (21)  and  (22))  contain  the  twelve  un¬ 
known  perturbation  quantities  p^,  pi,  u^,  vj,  or l,  Hi,  J3i,  jy,  jx,  ex, 
and  ey.  The  first  four  of  these  are  taken  as  the  primary  dependent 

variables  along  with  jx  and  jy.  In  order  to  obtain  a  closed  set  of  equa¬ 
tions,  additional  relations  are  needed  between  these  and  the  remaining 
perturbation  quantities. 

The  perturbed  state  properties  are  assumed  to  follow  an  ideal  gas 
behavior,  both  calorically  and  thermally.  Based  on  this  assumption, 
the  perturbed  enthalpy,  pressure,  density,  and  temperature  can  be 
related  in  the  following  manner: 

Y  PQ  Pi  Pi 

Hi  - - -  ( - __)  +  u  ui  (23) 

Y  -  IP  Po  *0  ° 

*  o 

and 


Ti  Px  Pi  (24) 

T“  "  IT  ”  p~ 

o  *o  Ho 

where  7  is  the  specific  heat  ratio  which  is  not  perturbed.  When  making 

calculations,  the  specific \he at  ratio  was  evaluated  at  the  thermodynamic 
conditions  given  by  the  steady-state  solution. 


The  electrical  conductivity,  <ji,  and  Hall  coefficient,  jBj,  will  now 
be  related  to  the  thermodynamic  parameters,  Tl,  p\,  pj,  Te^,  and 
nel.  The  Hall  coefficient  is  defined  as 


P  =  (JUT 


(25) 


where  e  is  the  electronic  charge,  me  is  the  electron  mass,  and  t  is 
the  mean  free-time  between  collision  of  electrons  with  other  particles. 
In  like  manner  the  electrical  conductivity  can  be  written  as 
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where  ne  is  the  electron  number  density.  Eliminating  t  between 
Eqs,  (25)  and  (26)  and  perturbing,  the  following  relation  between  cri, 
fil,  and  n^  results 


0i  01  nej, 

8  cr  n„ 
po  o  e0 


(27) 


Since  the  exact  expression  for  the  electrical  conductivity  is  more 
complex  than  is  necessary  for  this  analysis,  a  simpler  form  is  used 


a 


e 


Bi 

Te 


(28) 


where  Sm  is  the  seed  fraction  on  a  mole  basis  and  Te  is  the  electron 
temperature.  The  constants  a,  b,  and  Bj  can  be  picked  from  curve 
fits  of  the  electrical  conductivity  for  the  range  of  seed  fraction  and 
temperature  of  interest.  The  values  used  in  this  analysis  were 
a  =  0.  27,  b  =  0.  75,  and  Bj  =  22,  180°K.  These  correspond  to  a  seeded 
air  plasma  dominated  by  electron-neutral  collisions.  The  perturbed 
conductivity  can  now  be  written  as 


CTl  P1  1  Pi  Bi  V 

- a - i  —  +  (b  + 

ao  Pq  2  Po  Teo 


ex 


(29) 


where  it  has  been  assumed  that  the  perturbation  of  the  seed  fraction 
follows  the  gas  density;  i.e.,  (S^/S^)  =  ( Pl/pQ )■ 


The  electron  number  density  can  be  obtained  from  the  Saha  equa¬ 
tion  evaluated  at  the  electron  temperature 


n  2  r  2tt  in 

e  _  _ e 

n„  -  n  L  h2 

K.  e 

where  nj^  is  the  number  density  of  seed  nuclei  present  (seed  rate), 
k  is  Boltzmann's  constant,  h  is  Plank's  constant,  and  Vj  is  the  ioni¬ 
zation  energy  of  the  seed.  For  disturbances  of  reasonably  low  fre¬ 
quency  (up  to  about  a  megahertz),  the  electrons  are  assumed  to  re¬ 
main  in  local  Saha  equilibrium  and  from  Eq,  (30)  the  perturbed  elec¬ 
tron  number  density  can  be  expressed  as 


T  a/2  e" 
e 


e  V± 
Ft! 


(30) 


a. 


„  W3  6  Vi 
a  )  (—  + - 

°  2  k  T0 
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where  a0 , is  the  steady -state  ionization  fraction.  Using  Eqs.  (29) 
and  (31)  in  Eq.  (27)  results  in  the  following  expression  for  the  per¬ 
turbed  Hall  coefficient 


1  Pi 


- +  (a  - 

2  Po 


(32) 


’mo 


2  S, 


m. 


+ 


The  only  variables  remaining  unspecified  in  Eqs.  (29)  and  (32)  are 
the  steady-state  and  perturbed  electron  temperatures  (T6o  and  Tei). 

Since  the  purpose  of  this  analysis  is  to  investigate  the  effects  of  various 
parameters  on  the  stability  of  disturbances  in  the  flow,  the  steady- state 
electron  temperature  will  be  specified  by  fixing  the  quantity  (Teo  -  T0). 

When  the  electrons  are  in  thermal  equilibrium  (T0  =  TgQ)  Eq.  (24)  can 

be  used  for  the  perturbed  electron  temperature.  When  thermal  non- 
equilibrium  exists,  the  perturbed  electron  temperature  is  determined 
from  the  electron  energy  equation.  Simplifying  the  electron  energy 
equation  by  equating  the  energy  input  to  the  electron  gas  with  the  energy 
loss  attributable  to  collisions  of  electrons  with  other  species,  we  have 


where  6 
and  ves 


J2  J  . 

—  =  —  k  n,  m 

CT  2  e  e 


6  v 

Y  (  .  — )  <T„  -  T) 

s^e 


m_ 


(33) 


s 


s  is  the  energy  loss  factor  for  electron  collisions  with  species  s 
is  the  collision  frequency  of  the  electrons  with  species  s. 


An  effective  electron  energy  loss  factor  can  be  defined  as 


m. 


'eff 


r  6  v 

s4e  s  ms  es 

E  v 
s^e  es 


ms 

=  T  E  6_  —  v 

s  me  es 


s^e 


Then  using  Eq.  (26),  Eq.  (33)  is  written  as 

-  I  *  «eff(Te  '  T) 


(34) 
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Assuming  6eff  is  constant  and  using  Eqs.  (29)  and  (31)  the  perturbed 
equation  is 


2  J  j  +  2  J  j 
- O - -  (2  m  + 


ei 


Pi 


T  2  .2 
J  +  J 

X  V 


-  T 


-  T  p 
o  o 


(35) 


2  Sm 

+  (-  ° 


2  S  -  a.  T 
m©  o  e 


^ }  Pi 

- 


where 

,  S«*o  -  ao  w3  eVi 

|i  i  =  ( - )  (—  +  - ) 

2  S  -  a  2  kTe 
o  eo 

Equation  (35)  is  an  implicit  expression  for  Te^/TeQ  since  jy 

and  jx  depend  on  Tel/Teo  through  a\!a0  and  |3i//30.  An  explicit 

expression  will  be  obtained  for  the  various  cases  considered  in  the 
following  sections. 


The  number  of  unknowns  has  been  reduced  to  six:  pi,  pi,  u].,  vi, 
ey,  and  ex.  Additional  conditions  for  the  electrical  properties  must  be 
obtained  to  form  a  closed  set  of  equations.  In  addition  to  Eqs.  (17), 
(18),  (19),  and  (20),  the  perturbation  quantities  must  also  satisfy,  in 
some  sense,  the  following  relations 


-  dJx  dj 

7.  j - -  +  - -  =  0 

dx  dy 


(36) 


Vxe  =  0  or 


Sx 


de 

dy 


x 


(37) 


E 


y 


$  j 
oJy 


+  *  iJ, 


(38) 


(39) 
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where  4>i  and  Ii  are  the  perturbed  tangent  of  the  equipotential  slant 
angle  and  total  current,  respectively.  Although  the  equipotential  sur¬ 
face  may  remain  fixed  at  the  -walls  because  of  the  presence  of  the 
slanted  electrodes,  the  local  equipotential  surface  within  the  gas  may 
change  if  Ex  and  Ey  are  not  perturbed  in  such  a  manner  as  to  remain 
in  the  same  ratio  (ey  ^  <t>0ex). 

In  order  to  investigate  the  stability  of  the  slant-wall  MHD  acceler¬ 
ator  to  small  disturbances,  three  essentially  different  cases  will  be 
considered.  In  two  of  the  cases,  it  is  assumed  that  the  disturbances 
are  planar  waves  moving  in  the  axial  direction.  The  third  case  con¬ 
siders  disturbances  which  can  be  a  function  of  the  transverse  coordi¬ 
nate  as  well  as  the  axial  coordinate. 


SECTION  III 

ONE-DIMENSIONAL  PERTURBATIONS 


The  perturbations  are  assumed  to  be  in  the  form  of  one -dimensional 
plane  traveling  waves  moving  along  the  x  axis  and  are  dependent  only  on 
the  x  coordinate.  The  time  and  spatial  derivatives  of  the  small  disturb¬ 
ances  are  reduced  by  Fourier  analyzing  in  time  and  space  by  assuming 
the  waves  are  of  the  form 

ei(kx-ujt)  ^40) 

where  k  is  the  wave  number  which  is  assumed  to  be  real  and 
u  =  (Ji  +  iod2  is  the  complex  frequency.  The  wavelength  X  =  2irlk  and 
the  frequency  u  are  taken  to  be  independent  of  x.  This  requires  that 
the  wavelength  must  be  short  compared  to  the  distance  over  which  the 
steady-state  parameters  vary  appreciably. 

The  stability  of  the  disturbance  can  also  be  investigated  by  assum¬ 
ing  k  complex  and  id  real.  This  is  equivalent  to  investigating  the  growth 
of  the  disturbance  as  it  propagates  in  the  device  at  a  given  frequency,  in 
contrast  to  investigating  the  growth  in  time  for  a  given  wavelength  (k). 

It  was  shown  in  Ref.  10  that  a  wave  must  grow  in  time  and  space  simul¬ 
taneously  if  the  instability  is  indicated  in  either  case.  Equivalent  re¬ 
sults  were  found  in  the  cases  considered  here  using  either  approach. 

In  the  following,  k  is  assumed  real  and  id  is  assumed  complex. 
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When  the  perturbations  are  independent  of  y,  Eqs.  (36)  and  (40)  re¬ 
duce  to 


aj 


X 


dx 


0  or  j  =  0 
dx 


de 

— —  =  0  or  e„  =  0 
5x  ^ 


(41) 


This  implies  then  that  the  perturbation  current  flows  transverse  to  the 
wave  front,  and  the  electric  field  is  perturbed  normal  to  the  wave 
front.  The  plasma  is  infinite  compared  to  the  scale  of  the  disturbance, 
andEq.  (41)  represents  the  additional  relationship  needed  to  complete 
the  solution.  Equations  (38)  and  (3  9)  may  be  used  to  determine  the  local 
perturbations  and  Ii  but  are  not  necessary  for  completing  the  analysis. 
For  the  condition  in  Eq.  (41),  Eqs.  (38)  and  (39)  may  be  expressed  as 


#  i  ex 

*o  Ex 

ii  =  *ojy  +  ft  j.Jy 


(42) 


i.  e. ,  both  the  equipotential  slant  angle  and  total  current  are  changed 
from  their  steady-state  value. 

If  it  is  assumed  that  the  equipotential  angle  6  is  not  perturbed  and 
that  the  total  current  is  constant,  then  from  Eqs.  (38)  and  (39)  it  follows 
that 

ey  ~  ®oex  (43) 

Jx  -  -  *ojy  (44) 

This  implies  that  the  plasma  is  finite  and  the  accelerator  conducting 
walls  have  a  considerable  influence  on  the  disturbance.  These  relations 
together  with  the  assumed  form  of  the  wave,  Eq.  (40),  are  incompatible 
with  Eqs.  (36)  and  (37)  which  cannot  be  satisfied  exactly.  This  is  anal¬ 
ogous  to  the  steady-state  quasi-one-dimensional  solution  where  the 
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equations  are  not  satisfied  in  an  exact  sense.  The  essential  require¬ 
ment  is  that  J  and  E  are  slowly  varying  functions  of  x  so  that  the  one- 
dimensional  equations  are  approximately  satisfied.  The  same  require¬ 
ment  for  the  perturbation  quantities  ey  and  jx  approximately  satisfies 
the  perturbation  equations.  From  Eqs.  (37)  and  (40)  it  can  be  seen 
that 

de 

— -  =  i  k  e  «  0 

2>x  y 

implies  that  the  disturbance  must  necessarily  have  along  wavelength 
(k  <  1)  and/ or  the  slant  angle  6  is  small  so  that  ey  =  4>0  ex  is  small. 

The  requirement  that  the  disturbance  have  a  relatively  long  wave¬ 
length  contradicts  the  earlier  requirement  that  the  wavelength  should 
be  short  compared  to  the  x  variation  of  the  various  steady-state  param¬ 
eters  (uD,  pQ,  pQ).  It  can  be  seen  that  these  conditions  cannot  always 
be  satisfied  simultaneously  and,  hence,  all  of  the  equations  are  not 
satisfied  exactly.  However,  the  analysis  should  indicate  what  effect 
boundaries  have  on  disturbance  growth.  Equations  (43)  and  (44)  repre¬ 
sent  the  additional  relationships  needed  to  complete  the  solution. 

Hence,  two  different  cases  can  be  considered:  (1)  the  plasma  is 
infinite  compared  to  the  scale  of  the  disturbance,  and  Eqs.  (38)  and 
(39)  are  redundant;  (2)  the  plasma  is  finite  compared  to  the  scale  of 
disturbance,  and  the  boundary  conditions  of  constant  current  and  con¬ 
stant  equipotential  angle  must  be  satisfied. 


3.1  DISTURBANCE  IN  AN  INFINITE  PLASMA  (CASE  I) 

Since  the  perturbation  velocity  field  is  approximately  irrotational 
(vxv  =  0),  Ref.  8,  the  transverse  perturbation  velocity,  vi,  is  approxi¬ 
mately  zero.  Using  this  fact  and  the  results  from  Eq.  (41)  in  Eq.  (19) 
shows  that  the  perturbation  pressure  also  depends  only  on  the  x  coordi¬ 
nate.  The  y-rhomentum  equation  has  thus  been  approximately  satisfied. 
The  current  density,  jy,  joule  heating,  jyEy  +  Jxex,  conductivity,  cri. 
Hall  coefficient,  01,  and  electron  temperature  can  now  be  expressed  in 
terms  of  pi,  p and  uj  which  are  taken  to  be  the  independent  variables. 

Using  Eqs.  (21)  and  (22),  noting  jx  =  ey  =  0,  and  eliminating  ex 
yields  the  following  expression  for  jy 

ui  oi  8 1 

j  =  -  qa  —  +  q3  —  +  q4  a"  (45) 

J  o  o  wo 
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where 


q2  =  a  u  B 
o  o 


Qg  —  J  —  3  J 

y  OX 

q4  =  f3  J 
ox 


The  joule  heating,  <Jxex  +  jyEy,  can  be  expressed  as  follows 


Pi 


Jxex  +  Vy  -  rx  —  + 


Ui 


U 


(46) 


where 


rx  -  <Jy  -  soJx)Ey  -  jx*/aef 

e  j 

r*  «  -j-2  (Jy  +  eQJx  +  ®0e  ) 

o  * 


r3  =  -  B(a  E  +  M  ) 
o  o  y  Ko  x 


The  perturbed  quantities  ai/aQ  and  J3i/j30  can  now  be  related  to  ‘ 
the  variables  ui,  pi,  and  pi  by  using  the  expressions  developed  in  the 
preceding  sections.  Arriving  at  the  final  expressions  involves  some 
tedious  algebra,  and  their  development  can  be  found  in  Appendix  I 
along  with  the  final  coefficients.  Using  the  coefficients  from  Appen¬ 
dix  I,  jy  and  the  joule  heating  can  be  expressed  as 


3 


E 

y  y 


+  e  J 
x  x 


Pi 

Pi 

Ul 

—  + 

32 

—  + 

a3 

po 

Po 

uo 

Pi 

Pi 

Ui 

—  + 

b2 

—  + 

b3 

— 

Po 

Po 

uo 

(47) 


(48) 


The  perturbed  Eqs.  (17),  (18),  and  (20)  can  now  be  reduced  to 
three  homogeneous  algebraic  equations  in  the  unknowns  uj,  pj_,  and 
Pj.  Substituting  Eqs.  (47)  and  (48)  into  Eqs.  (17),  (18),  and  (20)  and 
using  Eqs.  (23)  and  (40),  the  algebraic  equations  can  be  written  in 
nondimensional  form  as 
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(uu .  -  k  +  iGl)  +  —  (-  k  +  iGi)  =  0 

uo 


Pi  Ulr  _  _  ,  P 1  k 

—  Ai  +  —  lA2  +  i(-  ID  +  k)  J+  —  (-  N2  +  i  = 

Pc  uQ  PQ  rM2 


—  [Bi  +  i(ou  -  k)]  +  — -  [B2  +  i (uu  -  k)  (y  -  1)M2] 
P  o  uo 

Pi  . 

+  —  [Ba  +  i( - +  k)]  =  0 

Po  Y 


where 


L  dpQ 
=  Po  dX 


G  2 


G3 


L_  duo 
uG  dx 

L_  d£o 

Po  dx 


A 1  =  G2  -  N 1 


A  2  —  G2  + '  N 


N  = 


Ni  = 


a  LB2 
o 

P  u 
Ko  o 

a  iBL 


N  2  = 


PqV 


(49) 

0  (50) 


(51) 
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Bi  =  —  G3  +  2Gi  +  ■■■'■  ■■  —  "  (J*  E  —  bi) 


y  0  u 
'  Ko  o 


B2  =  G2M2(r  -  1)  + 


r  ^  1  L 

r  p0u0 


(J •  E  -  b2) 


B3  =  -  G  i-l— 


-  1  L 


b3 


Y  Pnu 


00 


U)  =  UUl  +  1U)2  = 


luL 


k  -  kL 


and  L  is  a  characteristic  length  which  is  taken  to  be  one  meter  for  the 
cases  considered.  In  order  for  a  nontrivial  solution  to  exist  for  the 
set  (49-51),  the  determinant  of  coefficients  must  equal  zero.  This  re¬ 
sults  in  a  cubic  equation  in  the  complex  frequency  u,  i.  e. 


[■ 


i3  +  (-  3k  +  iCi)aua  +  k2(3 - )  +  C2  +  ikC3 

Ms 

+  P  (—  -  1)  +  +  ik2CB  +  ICb  =  0 

Ma 


]■ 


(52) 


where  the  C's  are  real  and  depend  on  the  steady- state  variables.  For 
a  given  steady -state  solution  satisfying  Eqs.  (6)  through  (8),  and  a 
given  value  of  k,  the  roots  are  determined  using  a  computer  program 
developed  at  AEDC  for  polynomials  with  complex  coefficients.  The 
criterion  for  instability  is  that  the  imaginary  part  of  the  root,  u,  be 
greater  than  zero;  i.  e. ,  W2  >  0*  This  implies  that  the  disturbance 
can  grow  with  respect  to  time  at  a  given  x  station  since 

i(kx-ctfit)  i(kx-u)it)  u)2t 
e  =  e  e 

It  can  be  shown  that  when  the  steady-state  gradient  terms  are 
zero  and  when  there  are  no  electrical  effects  present  in  the  perturbed 
equations,  there  are  three  real  roots  (Ref.  13)  of  Eq.  (52). 
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(u  -  a  )k 

0  0 

(53) 

u  k 

0 

(54) 

(u  +  a  )k 

0  0 

(55) 

The  root  (uQ  -  aQ)k  corresponds  to  a  wave  moving  upstream  relative  to 
the  gas  at  the  sound  speed  while  (uQ  +  aQ)k  is  a  wave  moving  downsteam 
relative  to  the  gas  at  the  sound  speed.  The  root  Ujjk  corresponds  to  a 
wave  moving  with  the  gas  and  is  identified  as  a  perturbation  in  enthalpy. 
Landau  and  Lifshitz,  Ref.  14,  show  that  the  entropy,  vorticity,  and 
enthalpy  perturbations  move  with  the  gas  velocity.  In  the  presence  of 
electromagnetic  effects,  the  disturbances  are  found  to  be  propagated  at 
the  speeds  implied  by  Eqs.  (53)  through  (55)  when  the  wavelength  is 
sufficiently  short.  These  will  be  referred  to  as  the  upstream,  down¬ 
stream,  and  enthalpy  modes.  The  disturbance  velocity  deviates  from 
these  values  as  the  wavelength  increases.  It  should  be  noted  that  the 
analysis  is  not  valid  at  extremely  short  wavelengths  because  of  the 
upper  limit  placed  on  the  disturbance  frequency  by  the  Saha  equilibrium 
assumption. 

Using  the  steady-state  solution  shown  in  Fig.  3,  the  dispersion 
relation,  Eq.  (52),  was  solved  at  various  points  within  the  channel. 

The  variation  of  u>2  with  x  is  shown  in  Fig.  4  for  the  enthalpy,  upstream 
and  downstream  wave.  Each  point  represents  the  local  value  of  the 
growth  rate  and  each  point  is  calculated  independent  of  the  growth  rate 
at  any  other  x  station.  It  can  be  seen  that  the  enthalpy  wave  exhibits 
the  largest  instability  (t02  >  1)  and  u>2  for  all  three  waves  shows  a 
marked  dependence  on  the  x  position  in  the  accelerator.  Each  x  station 
corresponds  to  a  different  orientation  of  the  wave  front  relative  to  the 
current  vector  as  well  as  different  values  of  the  steady-state  variables. 
The  dependence  of  wave  growth  on  wave  orientation  was  indicated  in 
Ref.  3  for  the  magneto -acoustic -type  wave.  It  can  be  shown  (Section  3.  3) 
that  the  maximum  joule  heating  occurs  when  the  angle  between  the  current 

vector  and  wave  front  is  given  approximately  by  e  max  =  ^  tan-1  |3.  A 

plot  of  the  angle  between  the  wave  front  and  current  vector  is  shown  in 
Fig.  5.  The  disturbance  is  oriented  near  the  angle  for  maximum  ampli¬ 
fication  at  the  channel  entrance  but  the  orientation  becomes  unfavorable 
for  growth  at  subsequent  downstream  stations.  The  decreasing  growth 
rate  (£2)  of  the  enthalpy  wave  with  increasing  x  (Fig.  4)  lends  credence 
to  the  dependence  on  orientation.  The  strong  variation  of  w  with  x  is 
apparently  attributable  to  the  variation  of  the  electrical  quantities 
rather  than  .the  x  dependence  of  the  gas  dynamic  or  thermodynamic  vari¬ 
ables. 
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x,  meters 


Fig.  4  Local  Disturbance  Growth  Rate  within  Accelerator  for  Case  I 
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MAXIMUM 


Fig.  5  Angle  between  Wave  Front  and  Current  Vector,  e  =  tan-1  (-  J„/Jy) 


21 


AEDC-TR-70-290 


In  order  to  determine  the  effects  of  the  steady- state  gradients,  the 
quantities  G^,  G2,  and  G3  were  set  to  zero  and  Eq.  (52)  was  solved  for 
this  case.  Neglecting  the  gradient  terms  increases  the  value  of  W2 
slightly  for  the  enthalpy  wave  but  does  not  change  the  general  character 
of  the  x  dependence  as  shown  in  Fig.  4.  Reference  13  found  a  decrease 
in  the  growth  when  gradients  were  neglected  in  the  generator  analysis. 
The  increase  in  growth  attributable  to  the  presence  of  gradients  was 
attributed  to  the  mechanism  described  by  Morse  (Ref.  15)  for  a  three- 
fluid  partially  ionized  plasma.  This  mechanism  depends  on  the  pres¬ 
ence  of  an  electric  field  either  parallel  or  antiparallel  to  a  density 
gradient  and  also  on  the  mobility  and  diffusion  of  the  electrons  and  ions 
within  the  gas.  This  mechanism  is  not  dominant  in  accelerators  or 
generators,  but  rather  the  gradients  contribute  convectively  to  the 
growth  or  decay  of  the  waves.  The  contributions  of  the  steady-state 
gradients  to  the  wave  growth  rate  were  shown  in  Ref.  8  to  be  approxi¬ 
mately 

4 

where  the  upper  sign  refers  to  downstream  waves  and  the  lower  sign  to 
upstream  waves.  The  contribution  to  the  amplification  or  attenuation  of 
a  given  wave  depends  only  on  the  sign  of  the  gradient. 

The  wavelength  dependence  of  u>2  the  nondimensional  wave 

Ml  ^WcLVGv 

speed  (—  = - )  are  shown  in  Figs.  6  and  7.  It  is  noted  that  the 

k  uo 

waves  travel  at  their  nondispersive  speeds  for  wavelengths  less  than 
approximately  0.  1  m.  The  wave  speeds  show  a  marked  deviation  from 
their  small  interaction  values  at  the  longer  wavelengths;  however,  at 
the  longer  wavelengths,  the  linearized  theory  breaks  down  and  the 
complex  frequency  w  can  no  longer  be  assumed  to  be,  even  approxi¬ 
mately,  independent  of  the  x  coordinate.  The  enthalpy  wave  u>2  is  noted 
to  be  almost  independent  of  wavelength  for  short  wavelengths,  but  then 
decreases  for  longer  wavelengths.  This  is  contrasted  to  the  results  of 
Ref.  6  where,  because  of  heat  conduction  losses,  the  growth  increases 
with  wavelength. 

The  effect  of  thermal  nonequilibrium  on  disturbance  stability  is 
investigated  using  the  perturbation  form  of  Te  given  in  Appendix  I. 

The  results  for  values  of  Te  up  to  200°K  higher  than  the  gas  tempera¬ 
ture  are  shown  in  Fig.  8.  The  nonequilibrium  effects  are  quite  large 
at  the  entrance  to  the  channel  where  the  enthalpy  wave  shows  the 
largest  value  of  £>2.  It  was  found  that  4>i  is  zero  for  Te  -  T  between 


1  +  (y  -  1  +  M0)  1  L  dpo  _  Mq  L  du{ 
1  ±  Mc  J  p  dx  2  u  dx 
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Fig.  6  Growth  Rate  versus  Wavelength  for  Case  I 


Fig.  7  Wave  Speed  versus  Wavelength  for  Case  I 
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Fig.  8  Growth  Rate  versus  Electron  Temperature  Elevation  for  Case  I 
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185  and  190°K  at  x  =  0.  02.  In  this  neighborhood  the  coefficients  of  the 
perturbed  electron  energy  equation  become  unbounded  and  the  linear¬ 
ized  analysis  breaks  down.  The  nonequilibrium  effects  are  quite  differ¬ 
ent  at  different  points  within  the  channel.  Reference  2  also  showed  that 
thermal  nonequilibrium  can  markedly  change  the  value  of  U2  f°r  the  up¬ 
stream  and  downstream  waves. 


3.2  DISTURBANCE  IN  A  BOUNDED  PLASMA  (CASE  II) 

The  stability  analysis  of  MHD  slant -wall -type  accelerator  flows, 
assuming  that  ey  =  jx  =  0,  indicates  that  the  enthalpy,  wave  can  become 
unstable  and  the  disturbance  may  become  large.  The  growth  rate  of  the 
enthalpy  wave  is  very  dependent  on  the  orientation  of  the  current  vector 
to  the  wave  front.  However,  different  assumptions  within  the  stability 
analysis  can  lead  to  radically  different  results.  If  different  boundary 
conditions  are  imposed  on  the  perturbed  electrical  quantities,  even 
though  the  same  steady-state  flow  exists,  different  stability  criteria 
will  result. 


Considering  the  case  where  the  plasma  is  finite  compared  to  the 
scale  of  the  disturbance,  the  current  density,  jy,  and  the  joule  heating, 
E0;j  +  e*  J0,  can  again  be  related  to  ui,  pi,  and  pi.  Making  use  of  the 
conditions  in  Eqs.  (43)  and  (44)  the  coefficients  for  the  jy  expression 
given  in  Eq.  (45)  are 


1  +  l0 2 

q3  =  -  q2 


q*  = 


1  +  *o  A 


The  coefficients  of  the  perturbed  joule  heating  term,  Eq.  (46),  become 

u  B(0  -  *  ) 


Yi 


T  r  u  tup  -  « 

. .  i  r.  +  — ° — °_i 
Aix  i  +  «*  J 


rs  = 


= -  r — - —  <|  -~2- .  p  u  b)| 

a  Ll  +  *  2  A  a0  o  o  J 
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I  uoB^o  -  V 
r3  A  x  +  i  2 


(56) 


The  coefficients  for  the  perturbed  electron  temperature  given  in 
Appendix  I  are  the  same  except  for  qi  which  becomes 


Qi  = 


2»y  -  »oV 
J  2  -  J  2 


The  development  of  the  dispersion  relation  follows  as  before  with  the 
dispersion  relation,  Eq.  (52),  as  the  result. 


Using  the  steady- state  solution  shown  in  Fig.  3,  tj  is  determined 
as  before  with  the  results  shown  in  Figs.  9  and  10.  It  is  shown  that 
the  enthalpy  and  downstream  waves  are  always  damped  whereas  the 
upstream  wave  has  a  positive  ujj.  The  corresponding  wave  speeds 
have  their  nondispersive  values  at  short  wavelengths  but  diverge  from 
these  at  longer  wavelengths.  The  effects  of  thermal  nonequilibrium 
are  shown  in  Fig.  11.  The  higher  electron  temperature  appears  to 
move  U2  toward  the  neutral  stability  point  (w2  =  0)  in  the  case  con¬ 
sidered. 


It  is  shown  that  instabilities  may  exist  with  the  wave  amplification 
depending  on  the  assumptions  used  in  the  analysis.  When  the  pertur¬ 
bations  in  Ey  and  Jx  were  zero,  the  enthalpy  wave  was  amplified  and 
would  grow  by  a  factor  of  105  in  one  meter  if  W2  remained  constant  at 
the  initial  value  of  12.  When  the  perturbations  in  Ey  and  Jx  were 
non-zero,  the  upstream  wave  was  amplified  and  would  grow  by  a  factor 
of  103  in  a  distance  of  one  meter  if  ^2  had  a  constant  value  of  six. 

These  growth  factors  can  be  taken  only  as  an  indication  of  the  possible 
amplification,  even  at  the  short  wavelengths,  since  contributions  from 
nonlinear  effects  could  become  large  and  «2  is  not  independent  of  posi¬ 
tion  in  the  channel.  Electron  thermal  nonequilibrium  is  shown  to  have 
a  marked  effect  on  the  growth  when  the  plasma  is  infinite  and  very 
little  effect  when  the  accelerator  walls  influence  the  disturbance.  Neg¬ 
lecting  contributions  from  the  steady- state  gradients  changed  the  mag¬ 
nitude  of  U2  only  slightly. 


3.3  APPROXIMATE  ANALYSIS 

Some  insight  into  the  contributing  factors  to  instability  can  be  ob¬ 
tained  from  an  approximate  analysis  similar  to  that  used  in  Ref.  3. 
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Fig.  9  Local  Disturbance  Growth  Rate  within  the  Accelerator  for  Case  II 
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Fig,  10  Growth  Rate  versus  Wavelength  for  Case  II 
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Neglecting  the  steady-state  gradient  terms,  the  perturbed  equations 
can  be  written  as 
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(59) 


where  a  different  form  of  the  energy  equation  is  used.  In  order  to  ob¬ 
tain  an  approximate  expression  for  id2,  the  right-hand  side  (RHS)  of 
Eqs.  (58)  and  (59)  are  computed  using  the  adiabatic  approximation 
Pi  Pi 

—  =  7  — .  Assuming  constant  collision  cross  sections,  the  perturbed 
Po  Po 

conductivity  and  Hall  coefficients  can  be  written  as 


Pi  r  +  i  p i 


oi  =  revi 

ac  L  2kT 


(y  -  l)  + 


y  -  3 


1- 

J  Po 


=  R 


Pi 


P. 


(60) 


Using  the  adiabatic  approximation  and  Eqs.  (57)  and  (60),  the  RHS  of 
Eq.  (58)  can  be  written  as 


BL  BL  f  uu  -  T 

— i  jv  = - ?  r  ^(“s— ) 

)^U  y  o  u  *  l  k 

o  o  o  o  '■ 


■)  +  q3RCT  - 


2  J  Po 


Pl  .  P  * 

F  — -  =  (Fr  +  iF  )  — 


(61) 
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The  RHS  of  Eq.  (59)  can  be  written  as 

(r- Dl  r  Po^o1 


a  p  _ 
oyo  o 


—  {  2F^~p-)[jy  +  (ii)j  l-  r  <j  *  +  j  2)T  — 
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o  r  ip 

Ko  Ko 


(62) 


Combining  Eqs.  (57)  through  (62),  the  condition  for  a  nontrivial  solu¬ 
tion  is 

r  +  (ip  -  .  o 

iu  -  k  k  k 


Approximating  (u  -  k)  in  the  term  iD/(u  -  £)  by  ±$c/M  (which  is 
regarded  as  an  iteration  in  solving  the  cubic  in  £3)  an  approximate 
expression  for  U2  can  be  written  as 


11)2 


^±“Fr 

2  y  2 


(63) 


The  upper  sign  corresponds  to  the  downstream  wave  and  the  lower 
sign  refers  to  the  upstream  wave.  Because  of  the  various  approxi¬ 
mations,  the  root  corresponding  to  the  enthalpy  wave  has  been  elimi¬ 
nated. 


The  expression  for  W2  when  ey  =  jx  =  0  can  be  written  as 
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(64) 


All  terms  in  Eq.  (64),  except  for  the  last  term,  are  identified  in 
Eqs.  (39)  and  (40)  of  Ref.  3.  The  last  term  is  a  coupling  term 
required  by  the  presence  of  two  components  of  current  in  the  joule 
heating  term,  Dp.  A  plot  of  Eq.  (64)  using  the  steady-state  solution 
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of  Fig.  3  is  shown  in  Fig.  12a.  Comparison  of  Fig.  12a  with  Fig.  4 
for  the  upstream  and  downstream  wave  shows  that  the  simplifying 
assumptions  and  approximations  made  in  arriving  at  Eq.  (64)  affect 
the  magnitude  of  U2  slightly  but  do  not  alter  the  character  of  the 
x  dependence.  The  growth  or  decay  of  the  waves  is  shown  to  be  very 
dependent  on  the  magnitude  (and  sign)  of  Jx  and  Jy.  The  largest  con¬ 
tribution  to  instability  comes  from  the  third  term  which  contains  the 
axial  current  density,  Jx.  As  in  previous  analyses  (e.  g. ,  Refs.  3 
and  7)  the  wave  propagating  antiparallel  to  the  current  vector  is  ampli¬ 
fied.  Since  Jx  goes  through  zero,  the  downstream  wave  is  amplified 
near  the  channel  entrance,  but  the  upstream  wave  is  amplified  near 
the  exit. 


The  angle  between  the  current  vector  and  the  wave  front  (see 
Fig.  5)  which  gives  the  maximum  growth  rate  can  be  calculated  from 
Eq.  (64).  The  body  force  terms  (second  and  third  terms  of  Eq.  (64)) 
are  a  maximum  when 


tan  e 


Ra  +  — g— 

e( - ) 

Ra  -  2(r-l) 


The  joule  heating  terms  (last  two  terms  of  Eq.  (64))  are  a  maximum 
when 


tan  2e  -  P  ( - -  ■  ) 

ft 

a 


The  maximum  contribution  from  the  body  force  terms  occurs  when  the 
angle  between  the  current  vector  and  wave  front  is  approximately 
86  deg.  Since  the  force  terms  are  larger  than  the  joule  heating  terms, 
this  orientation  will  give  the  maximum  growth  rate.  As  can  be  seen 
from  Fig.  5,  this  optimum  orientation  is  not  present,  and  axially 
moving  acoustic  waves  are  not  significantly  amplified. 


The  expression  for  U2  when  ey  =  <t>0ex  and  jx  =  -<t>ojy  can  be 
written  as 
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(65) 


The  first  three  terms  are  related  to  fluctuations  in  the  forcing  term, 

Fr,  and  the  remaining  terms  are  related  to  changes  in  the  joule  heat¬ 
ing.  A  plot  of  t02  versus  x  for  this  case  is  shown  in  Fig.  12b  using 
the  steady-state  solution  of  Fig.  3.  The  result  is  very  similar  to  that 
obtained  by  solving  the  cubic  dispersion  relation  (Fig.  9).  The  largest 
contribution  comes  from  the  perturbation  in  the  forcing  term,  Fr,  i.  e. , 
the  second  and  third  terms  of  Eq.  (65).  The  body  force  term  is  inde¬ 
pendent  of  the  orientation  of  the  wave  front  to  the  current  vector.  The 
contribution  from  joule  heating  decreases  as  x  increases  and  is  negli¬ 
gible  at  the  end  of  the  channel. 


The  growth  of  the  enthalpy  wave  can  be  calculated  approximately 
by  considering  the  growth  of  a  stationary  enthalpy  disturbance  attrib¬ 
utable  to  joule  heating.  It  is  assumed  that  the  disturbance  is  a  plane 
wave  so  that  ey  =  jx  =  0  corresponding  to  Case  I.  The  perturbed 
energy  equation  takes  the  simple  approximate  form 
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(66) 


Using  Eqs.  (9)  and  (10),  the  perturbed  current  density  can  be  written 
as 


a  1  a  t  \  ^ 1 

jy  =  — (Jy  “  +  —  ^o^x 

*o  ^o 


(67) 


It  is  assumed  that  the  electrical  conductivity  and  Hall  coefficient  can 
be  written  as  power  law  functions  of  the  enthalpy;  i.  e. 


Wl 


(68) 
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(69) 


Using  Eqs.  (66)  through  (69),  the  perturbation  energy  equation  can  be 
written  as 
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(70) 


using  perfect  gas  relations. 

The  right-hand  side  of  Eq.  (70)  has  been  nondimens ionalized  using 
the  free -stream  velocity  so  that  it  now  corresponds  to  the  temporal 
growth  rate,  £2.  It  is  shown  that  the  first  term  of  Eq.  (70)  is  always 
destabilizing  since  Jy  >  Jx,  and  the  second  term  is  destabilizing  when 
Jx  is  negative  since,  in  general,  wi  >  W2-  Using  wi  =  6  and  W2  =  1, 
the  approximate  enthalpy  growth  rate  is  calculated  from  Eq.  (70)  using 
the  steady -state  solution  of  Fig.  3.  The  growth  rate  is  shown  in 
Fig.  13.  A  comparison  of  Fig.  13  with  Fig.  4  shows  that  Eq.  (70) 
underestimates  the  value  of  the  growth  rate  at  the  accelerator  entrance 
but  does  qualitatively  predict  the  growth  rate  distribution  through  the 
accelerator. 

The  angle  for  maximum  enthalpy  growth  can  be  found  from  Eq.  (70) 
by  relating  the  current  components  to  the  angle  between  the  wave  front 
and  current  vector  (defined  in  Fig.  5).  The  current  components  can  be 
written  as 

Jy  =  J  cos  e  (71) 


J  =  -  J  sin  e 
x 


(72) 


Using  Eqs.  (71)  and  (72)  in  Eq.  (70),  the  approximate  expression  for 
the  growth  rate  can  be  written  as 

L 
uo 

(73) 


L  J2  y  -  1 
—  (— )(- - )  W! 


u  CT 

O 


r  p 


{ 


cos  2e  +  sin  2e(l 


-  — )  1 

J 


35 


A  EDC-TR -70-290 


Fig.  12  Approximate  Growth  Rate  of  Magneto-Acoustic  Waves  as  a  Function  of  x 
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The  angle  for  maximum  growth  is  then  found  from  Eq.  (73)  as 

tan  2a  =  S0<1  -  (74) 

where  the  term  W2/wj  represents  the  contribution  from  the  perturbed 
Hall  coefficient  which  is  stabilizing.  The  orientation  most  favorable 
for  growth  is  near  the  channel  entrance,  and  the  wave  orientation  be¬ 
comes  less  favorable  for  growth  as  x  increases.  This  trend  agrees 
with  the  growth  rates  predicted  from  the  Case  I  analysis. 

The  approximate  results  shown  here,  together  with  the  results 
from  the  previous  cases  (ey  =  jx  =  0  and  ey  =  ^Vjex;  jx  =  -  $>0jy)  indi¬ 
cate  that  instabilities  may  exist  in  slant-wall  accelerators.  Although 
the  linearized  analysis  is  inadequate  for  predicting  the  disturbance 
growth  and  is  restricted  to  short  wavelength,  insight  is  gained  into 
the  instability  mechanism.  It  is  shown  that  different  assumptions 
concerning  the  electrical  boundary  conditions  give  radically  different 
results  for  stability  criteria.  In  an  attempt  to  remove  some  of  the 
restrictions  in  the  analysis,  a  two-dimensional  perturbation  model  is 
considered  in  the  next  section. 


SECTION  IV 

TWO-DIMENSIONAL  ANALYSIS 


The  perturbation  variables  must  satisfy  Eqs.  (17)  through  (20) 
together  with  the  conditions  imposed  by  Eqs.  (36)  through  (39).  Pre¬ 
viously,  it  has  been  assumed  that  the  imposed  disturbances  are  axial 
traveling  waves  which  depend  only  on  x  and  t.  This  restriction  will 
be  relaxed  to  investigate  the  stability  of  traveling  waves  which  are 
dependent  on  the  transverse  coordinate,  y,  as  well  as  x  and  t.  This 
approach  permits  a  more  realistic  treatment  because  the  perturbed  y 
momentum  equation  is  not  neglected  and  the  waves  are  not  restricted 
to  propagate  only  in  the  x  direction. 

The  conditions  Vxe  =  0  and  V x  v  =  0  imply  that  the  perturbed 
electric  field  and  velocity  are  derivable  from  potential  functions. 
Hence 


~  =  -  7V 
-  ?qp 


(75) 
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It  is  assumed  that 


V  ~  e 
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(76) 


The  components  of  perturbed  electric  field  and  velocity  can  now  be 
expressed  as 


e  =  -  ik  V 
x  x 


ey  "  -  lkyV 


iu  =  ikxcp 


vi  =  ikycp 


(77) 


The  conservation  of  current,  v  j  =  0,  implies  that  j  can  be  expressed 
as 


dy 

_  dj, 

dx 


(78) 


McCune  (Ref.  3)  usedEqs.  (75)  and  (78),  together  with  Ohm's 
law  and  the  momentum  and  continuity  equations,  to  determine  a  set  of 
differential  equations  which  together  with  appropriate  boundary  condi¬ 
tions  determine  V,  <f> ,  ijj,  and  p.  The  energy  equation  was  eliminated 
by  noting  that  dissipation  terms  are  much  smaller  than  the  body  force 
terms  for  the  case  considered  and,  hence,  the  adiabatic  approximation 

(Zi  =  7  fli)  can  be  used.  The  adiabatic  approximation  is  not  made  in 
\Po  Po/ 

this  analysis  and  conservation  of  current  is  not  satisfied  exactly. 

Ohm's  law  is  used  to  establish  the  relation  between  $  ,  V,  and  p. 
Equations  (17)  through  (20)  can  in  principle  be  solved  for  <j> ,  V,  p, 
and  p  and  then  ^  can  be  determined  from  Ohm's  law. 


The  assumed  form  of  the  perturbations,  Eq.  (76),  implies  that  the 
medium  is  infinite,  compared  to  the  scale  of  the  disturbance,  in  both  x 
and  y  direction.  It  is  not  possible,  then,  to  treat  reflection  or  attenua¬ 
tion  of  a  disturbance  at  fhe  accelerator  walls.  However,  the  initial 
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development  of  the  disturbance  can  be  investigated  and  the  analysis 
can  be  used  to  determine  if  the  necessary  conditions  for  instability  are 
present. 

The  perturbed  Ohm's  law,  Eqs.  (21)  and  (22),  can  be  .written  as 
=  —  =  -  ik  a  J  +  ik  8  <t  „V  +  ik  a  -Bcp 


x 


Sy 


ef  o  x~ 


Vo~ef'  T  y  ef* 

x  a"  -  &ef'“y  T  Ko-X'p 


x  ef 

a,  3,  <79> 

+  Bk  <p  +  J  - - 


jy  -  -  -  -  iky°efV  -  lkxBael<l>  -  ikx^oCTefV 

+  ^ky^o^°ef^  +  Jy  “  +  pef(,x  '  Vy^  <8<>) 

o  ° 

The  expressions  for  jx  and  jy  are  used  together,  with  the  perturbed 
joule  heating 


(J«E)i  =  Je  +Je  +  E  j  +  E  j 
xx  yy  xx  yy 


(81) 


and  the  expressions  given  in  Appendix  I  for  the  perturbed  conductivity 
and  Hall  coefficient  to  express  Eqs.  (17)  through  (20)  in  terms  of  <t> , 
V,  p,.  and  p. 

The  perturbed  equations  can  now  be  written  as 
Pi  r  Pi 


^  f  -  I  P1  f  -  -  TV 

—  +  1  -  Aa  +  ik  PV  —  +  J  i(N  -k  +  k  M  J  )  — 

Po  l  x  J  P0  L  ef  y  x  o  ef  J  uQ 

+  J  k  m  -  k  2  -  i(R  N  _k  -  k  A  )  V  — —  =  0 

\  x  x  vpo  ef  y  x  u  l 


(82) 


^  +  +iV}^+  {1(PoNefky  -  kx"ef>} 

+  {“*y  -  kxky  +  i(kyNef  +  kxEa)}  ^ 


u  BL 
o 


(83) 
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k  2  -  k  2  +  ik 

A  J 


lkxGl}-  0 


{b>  +  i(.- !.>}£♦{*  ♦*»,">}£ 

+  {-  i(kyb*  +  ^xVjs-gL  +  {<V  -  Fx2)  M2(r-D 

-  -  'I  cpi 

-  i(kyb6  -  kxB2)  j—  =  O  (85) 


where  the  following  new  definitions  are  made 

CTef  B=L 


J  BL 

F  =  -1 _ 

y  Pouo2 


Jx  BL 

* 


"  G2  -  Fy§a  -  ^3Pef(Fx  -  PQFy) 


A2  =  Ga  +  Nef  -  Fy?a  -  IUP0f(Fx  ~  P/y) 


Aa  =  V1+tiaPef(Fx-  W 


E 1  =  ?aFx  “  ^3  ^ef (Fy  +  PoFx) 


Ea  =  ^oNef  +  V3  “  ^ef(Fy  “  W 


Es  “  ^Fx  -  Ha3ef(Fy  +  3oFx) 
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bi  = 


bp  = 


b3  = 


b4  = 


b=  = 


b6  = 


1=1 

L 

Y 

P  u 
^0  0 

r- 1 

ILL 

Y 

p0  t 

1=1 

L 

Y 

pouo 

Y- 1 

Llm 

Y 

po  t 

Y-l 

L 

Y 

P  u 

0  0 

Y-l 

Y- 

Po  1 

2G 1 

-  G3 

(5  2  -  0oPefH3)J-E  +  p_p|l3(ETJv  - 


ef^  V~y"x 


E  J  ) 
x  y 


} 


("1  -  +  Pef|l2(EyJX 


ExJy) 


} 


J  B  +  g  ^BCp  E  -  E  ) 
y  ef  o  x  y 


} 


I  „(E  J 
ef  y  x 


E  J  ) 
x  y 


+  “ef  B  ao(EoE*  -  Ey> 


} 


ef  B(Ex  +  0 


E 

o  y  J 


—  —  y-\  L 


r  p0uo 


-)  -  bi 


_  — »  — »  y- 1  L 

Bp  =  M2  (y-l)Gp  +  J°  E(- - )  -  bB 

Y  P0uo 


B3  =  -  G 1  -  b3 


(86) 


The  dispersion  relation  is  obtained  by  setting  the  determinant  of 
coefficients  of  the  algebraic  set  of  equations,  Eqs.  (82)  through  (85)  to 
zero.  After  considerable  algebraic  manipulations,  the  resulting  dis¬ 
persion  relation  is  of  the  form 

CAiuu3  +  (CA2  +  iCBp)uu2  +  (CA3  +  iCB3)uj  +  CA4  +  iCB4  =  0 

(87) 
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The  dispersion  relation  is  again  a  cubic  and  in  the  absence  of  elec¬ 
tromagnetic  effects  reduces  to  the  equations  for  traveling  acoustic  and 
enthalpy  waves.  The  steady-state  solution  shown  in  Fig.  3  is  utilized 
with  the  investigation  restricted  to  thermal  equilibrium;  that  is,  Te  =  T. 

It  is  seen  from  Fig.  14  that  for  wavelength  in  the  y  direction  greater 
than  0.  1  m  and  for  large  values  of  kx,  the  growth  rate  reduces  to  that 
obtained  from  the  one -dimensional  analysis  with  ey  =  ex  =  0  (Fig.  4). 
This  was  expected  since  as  ky-*-  0  the  gradients  in  the  y  direction  are 
eliminated.  Also,  in  both  cases  the  perturbation  of  slant  angle  and 
total  current  are  nonzero,  whereas  in  the  case  where  ey  =  <fc0ex  and 
jx  =  -  4>oJy>  they  were  assumed  zero.  However,  the  two-dimensional 

perturbation  analysis  has  0,  but  this  finite  perturbation  in  the  axial 
current  does  not  change  the  stability  results  for  short  wavelengths. 


The  growth  rate  for  the  various  waves  is  plotted  versus  axial  wave¬ 
length  in  Fig.  15  and  versus  ky  in  Fig.  16.  The  growth  rate  for  the 
enthalpy  wave  is  shown  to  decrease  for  increasing  Xx  and  decreasing  Xy. 
The  growth  rate  for  the  upstream  wave  increases  slightly  for  increasing 
Xx  and  decreasing  Xy,  but  the  opposite  is  true  for  the  downstream  wave. 
Changing  Xx  (or  Xy)  is  equivalent  to  changing  the  wave  orientation  in  the 
flow,  with  the  angle  between  the  current  vector  and  wave  front  given  by 


tan  e 


1  lx 

’  ‘  S  s’ 


Then,  at  the  accelerator  entrance,  increasing  Xx/Xy  corresponds  to  a 
decreasing  angle  and  a  less  favorable  condition  for  amplification  of  the 
enthalpy  wave. 


The  wave  speed  dependence  on  wavelength  is  shown  in  Fig.  17.  As 
before,  the  waves  deviate  from  their  nondispersive  velocities  at  wave¬ 
lengths  (Xx)  on  the  order  of  0.  1  m.  The  wave  speeds  of  the  upstream 
and  downstream  waves  deviate  at  smaller  values  of  Xx  as  the  wavelength 
in  the  y  direction  decreases  below  one  meter. 
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SECTION  V 
SUMMARY 


The  linearized  stability  analyses  indicate  that  axially  moving 
magneto -acoustic  and  enthalpy  waves  can  be  amplified  in  slant-wall 
accelerators.  The  stability  and  degree  of  amplification  of  any  wave 
depend  on  the  electrical  boundary  conditions  and  wave  form  utilized. 

The  imaginary  component  of  the  complex  frequency  (£2,  which  indicates 
the  "degree"  of  instability)  was  found  to  be  a  strong  function  of  the  x 
coordinate.  The  strong  x-dependence  was  seen  to  be, related  to  the 
orientation  of  the  wave  with  respect  to  the  current  vector.  The  marked 
x  dependence  of  U2  and  the  magnitude  of  the  gradients  of  p0,  uo,  pQ 
indicate  that  the  analyses  are  restricted  to  wavelengths  of  0.  1  m  or  less. 
The  disturbance  traveled  at  the  local  sound  speed  (relative  to  the  flow) 
or  at  the  local  fluid  velocity  for  short  wavelengths,  but  lost  this  char¬ 
acter  as  the  wavelength  increased.  In  some  cases  (at  wavelengths 
greater  than  0. 1  m)  the  upstream  wave  moved  at  a  velocity  greater 
than  the  fluid  velocity,  and  thus  would  appear  to  have  a  negative  veloc¬ 
ity  relative  to  a  fixed  observer.  However,  the  validity  of  the  analyses 
is  doubtful  at  the  longer  wavelengths. 

Table  I  (Appendix  II)  summarizes  the  results  of  the  linearized 

analyses.  The  quantity  ew2T  corresponds  to  the  growth  of  a  disturbance 
in  the  device  one  meter  in  length  with  u 2  assumed  constant.  .  This  only 
provides  an  indication  of  the  possible  growth  of  a  disturbance  since  non¬ 
linear  effects  will  become  important  when  the  disturbance  becomes 
large.  It  is  predicted,  from  both  the  one -dimensional  and  the  two- 
dimensional  analyses,  that  the  enthalpy  wave  will  have  the  largest 
growth. 

The  enthalpy  wave  is  more  prone  to  becoming  unstable  in  a  large 
device  (ey  =  jx  =  0  or  2-D  analyses)  where  the  boundaries  have  no  influ¬ 
ence  on  the  disturbance.  The  disturbance  oriented  at  the  most  favor¬ 
able  angle  for  amplification  would  grow.  This  is  not  unlike  the  growth 
of  ionization  waves  with  the  resulting  striations.  The  upstream  acoustic 
wave,  on  the  other  hand,  is  more  apt  to  be  amplified  in  small  devices 
where  the  boundaries 'have  a  marked  effect  on  the  disturbance  (ey  =  <&c>ex 
and  jx  =  -  ^ojy).  In  this  case,  the  upstream  wave  is  amplified  regard¬ 
less  of  orientation. 

The  linearized  anlaysis  indicates  that  instabilities  can  exist  in  MHD 
accelerators,  but  the  analysis  cannot  predict  the  effect  on  device  per¬ 
formance.  A  linearized  analysis  cannot  consider  the  long  wavelengths 
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and  aperiodic  waves  which  may  be  present  in  MHD  devices  (Refs.  16 
and  17).  A  more  complete  analysis  which  considers  the  perturbed 
differential  equations  subject  to  the  proper  boundary  conditions  is 
necessary  to  predict  the  effect  of  disturbance  growth  on  accelerator 
performance. 
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APPENDIX  I 

CURRENT  AND  JOULE  HEATING  PERTURBATIONS 


In  order  to  relate  the  perturbed  transport  quantities  (0i  and  ai)  to 
the  perturbed  thermodynamic  properties,  it  is  necessary  to  use  the 
perturbed  electron  energy  equation  (Eq.  (35)) 


2  J  j  +  2  J  j  T~  Tei 

—111 - „  <2|ll  + - 2°_)  — 


pi 


T  2  T  2 

J„  +  J 

x  y 


Teo  "  To  Teo 


-  T_  p. 


+  (- 


2  S, 


m. 


2  Smn  "  a. 


Te, 


-  T 


a-i) 


The  perturbed  Hall  coefficient  and  electrical  conductivity  are  related  to 
p i,  pi,  and  Te  by  the  following  expressions 


Bi  .  Tei 


ai  P1  1  PI  *  V 

—  =  a - - +  (b+  - ) 

a  pQ  2  PQ  T  Te 


(1-2) 


3  x  1  Px 
—  =  - - +  (a  - 


Px 

-)  —  + 


Po 


2  P, 


2  Smo"  ao  Po 


r  B- 
{b  +  ~ 


S_  -  an  q  dVi 
.  (  ;°)(-  +  — 
2  S  -  a  2  kTe 


\  Tei 
'/  57 


(1-3) 


Using  Eqs,  (1-1),  (1-2),  and  (1-3),  an  explicit  expression  for  Tei/Te(, 
can  be  obtained 


■ai  ui  pi  pi 

-  -  $p2  -  +  cp3  -  +  Cp4  - 

u_  0  _  D 


(1-4) 
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where 


cpi  =  - 
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Tp  -  T 
eo  o 


+  pi(2  +  q  xq4  ) 


Bi 

-  qi(q3  +  q4)(b  +  — )  r 

eo  J 


cp2  =  qiq2/cpi 


,3  -  if 

3  <p,  V 


Sm, 
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(2  +  q^)  -  aq^qs  +  q4)  J 


1  r  To  qx  1 
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2  J, 


qi 


j  2  +  J  * 
X  y 


;  q2  «=  cjoUqB;  q3  «  J  -  PQJX»  q4  *  30J, 


Having  determined  the  temperature  perturbation,  the  expressions 
for  the  conductivity  and  Hall  coefficient  can  now  be  written  as 
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Pi  Pi 
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where 


=  _  I  +  cp4(b 


?3 


42 


43 


=  a  +  cp3  ( b  +  = — ) 

Bi 

=  cp2(b  +  — ) 

T 

e0 

=  ^  1  -  4  i<p4 

=  ?2  -  4  i<p3  -  Smo/(2  Smo 


44  =  ^3  -  41CP4 

Substituting  Eqs.  (1-5)  and  (1-6)  into  Eqs.  (45)  and  (46)  gives  the  follow¬ 
ing  expressions  for  the  current  and  joule  heating: 


(1-7) 
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where 


b  1  =  Yi$  2  +  Y  243 
b2  =  yz  +  YiAs  +  Y  244 
b3  =  Y  i§  1  +  Y  24  2 


a  1  =  t  2q3  +  43  q4 
a  2  =  5  iq3  +  4  2Q4 

as  =  “  <12  +  ?3Q[3  +  44^4 
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APPENDIX  II 
TABLE  I 

SUMMARY  OF  RESULTS 


Analysis 

x,  m 

Enthalpy 

Upstream 

Downstream 

One-Dimensional 

0.0 

Amplified  e^2^  «  106 

Attenuated 

Neutral  u)2  »  0 

Perturbation 

Case  I 

0.76 

Neutral  uus  0 

Attenuated 

Attenuated 

e  =  j  =0 
ey  Jx 

2.4 

Attenuated 

Amplified  e^2^  «  25.0 

Attenuated 

One-Dimensional 

0.0 

Attenuated 

Amplified  eU)2T  w  40 

Attenuated 

Perturbation 

Case  II 

0.76 

Attenuated 

Amplified  ews  «  300 

Attenuated 

ey  =  *oex 

2.4 

Attenuated 

Amplified  e^2*1  «  300 

Attenuated 

"  “  #0Jy 

Two-Dimensional 

0.0 

Amplified  e^2*1  »  10s 

% 

Attenuated 

Amplified  e^2^  »  4.4 

Perturbation 

X  =  1  m 

0.76 

Amplified  u)a  0.016 

Attenuated 

Attenuated 

y 

2.4 

Attenuated 

Amplified  eu’2T  m  7 

Attenuated 

T  -  Time  spent  in  device  one  meter  in  length 
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